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Abstract 

We study the behaviour of the contact process with rapid stirring 
on the lattice Z"^ in dimensions d> 3. This process was studied earher 
by Konno and Katori, who proved results for the speed of convergence 
of the critical value as the rate of stirring approaches infinity. In this 
article we improve the results of Konno and Kattori and establish the 
sharp asymptotics of the critical value in dimensions d > 3. 
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1 Introduction and the main result 



In this paper, we are going to study the behaviour of the so-called contact 
process with rapid stirring (see [5], [6]). The process is also known as a con- 
tact process combined with an exclusion process (see [ZIIH])- 

DeMasi, Ferrari and Lebowitz (see |2]) studied interacting particle systems 
on a lattice under the combined influence of spin flip and simple exchange 
dynamics. They proved that when the stirrings occur on a fast time scale of 
order the macroscopic density, defined on spatial scale e~^, evolves ac- 
cording to an autonomous nonlinear diffusion-reaction equation. Using the 
connection between a convergent sequence of such particle systems to a solu- 
tion of a react ion- diffusion equation, found by DeMasi, Ferrari and Lebowitz, 
Durrett and Neuhauser (in |3]) proved results about the existence of phase 
transitions when the stirring rate is large that apply to many different sys- 
tems. 

The contact process with rapid stirring was studied by Konno (see [6]), who 
described it via a system of interacting particles, on a lattice Z"^ (while for 
the proofs the process on the rescaled lattice Z'^/a/ZV was considered). The 
state of the process at time t is given by a function t]^ : Z'^ — )■ {0, 1}, where 
the value of rjf{x) is determined by the number of particles present at x at 
time t. In this setting particles die at rate 1, and give birth, onto the closest 
neighbouring sites, at rate A. In addition, values of rj^ at two neigbouring 
sites are exchanged at rate (stirring rate), and all the above mechanisms 
are independent. The primary goal of [6] was to improve the result of Durrett 
and Neuhauser, who showed the following: 

Theorem 1.1 (Durrett, Neuhauser, jS]) Let {r]f}t>o be the set-valued 
contact process with stirring, with the dynamics as described above, starting 
with a single occupied site at the origin. Let Q^o = {f]? 7^ 0? Vt > 0}, 
= 'P [fioo] , CLnd let 

Ae(iV) = inf{A > : > 0}. 

As N ^ oo, 

a. A,(iV) ^ 1. 

b. // A > 1, then ^ (A - 1)/A. 
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Konno used the methods of [T] to get a more detailed picture of the critical 
value Xc{N) as the stirring rate approaches infinity. The main result of [B] 
can be stated as follows. 



Theorem 1.2 (Konno, [6]) For all x > 0, let 

( l/xV3, d=l, 

(fd{x) = < logx/x, d = 2, 
[ 1/x, d>3. 

Then we have 

X,{N) - 1 ^ C^MN) 

where ~ means that if is small (large) then the right-hand side of the 
above is a lower (upper) hound of the left-hand side for large enough N . 

Moreover, Theorem 11.21 was refined in dimensions > 3: with the help of 
result of Katori (see |5]), it was shown in [6] that 

< liminf iV(Ac(iV) - 1) (1.1) 

A''— >-oo 

< limsupA^(Ac(iV) - 1) < 



(2d)(2rf-l) - N^oo 

G(0,0) - 1 



where G(-, ■) is the Green's function for the simple random walk on Z . 

The main goal of this article is to show that, in fact, the lower bound 
in (11.11) can be improved to the value '^^^2d~^ ' this, we get the sharp 
asymptotics of the critical value Xc{N). Before we state our main result, we 
need another piece of notation. Let Srf(x) = {y G Z*^ : \\y — x\\i = 1} denote 
"the neighbourhood of x in Z"', excluding x" . Then = 9^(0) will denote 
the neighbourhood of the origin. 

The main result of this article is 
Theorem 1.3 Let d>3. Then 

Xc{N) - 1 ~ ^, as N ^oo, 
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where 

^ oo 

^ = i^Ep[^"^s^]' (1-2) 

n=l 

{Ki}n>o is a symmetric random walk on if" starting at the origin, and ~ 
means the ratio approaches 1, as N approaches oo. 

To connect the result in the theorem with fll.ip let us state a simple 

Lemma 1.4 

2di9 = G(0,0) - 1. 

Remark 1.5 The lemma implies that our sharp asymptotics for the critical 
value coincides with the upper bound for the critical value in U.l\) . 

Proof of Lemma 11.41 By the Markov property 

^ oo 

G(0,0) = 1 + _^P[K-iG5,] 

n=l 

oo 

n=l 

where the last line follows since Vq = 0. □ 

Let us say a few words about the proofs. The structure of Konno's proofs 
of Theorem 11.21 follows the ideas of Bramson et alii (see pQ), who studied the 
long range contact process (LRCP) in a limiting regime, when the range M 
of the contact process, goes to infinity. The set up of Konno is very similar 
to that of Bramson et alii in [1], with the difference that the stirring speed 
goes to infinity, and not the range. Note that in pQ, the authors were able to 
prove an asymptotic result for LRCP which was later improved by Durrett 
and Perkins (see j^) where a sharp asymptotics for the convergence of the 
critical value was obtained for dimensions d >2. To prove the sharp asymp- 
totics, Durrett and Perkins, in [1], had further rescaled space and time and 
proved week convergence of the rescaled processes to super-Brownian motion 
with drift. This convergence almost immediately gives the lower bound for 
the critical value for LRCP. As for the upper bound, Durrett and Perkins 
bounded the LRCP from below by an oriented percolation process (for this 
they used again convergence to super-Brownian motion), and, by this, the 
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upper bound for the critical value was derived. 



As for the proof of our main result — Theorem 1 1.31 — it follows immediately 
that from (11.11) and Lemma [1.41 that it is sufficient to prove just the lower 
bound for the critical value. This makes the proofs far less complicated than 
those in In fact, we prove our result without proving a weak convergence 
result of rescaled processes to super-Brownian motion, which was one of the 
main technical ingredients of the proofs in 

The rest of the paper is organized as follows. Formal definitions of the 
contact process and "speeded- up" contact process are given in Section [2j 
Theorem 11.31 is proved in Section [21 

2 Formal definitions 

Before we proceed to the proofs, let us give formal definitions in this section. 

The contact process with rapid stirring takes place on the lattice . Fix 
parameter Q for this process. The state of the process at time t is given 
by a function j]^ : if- — )■ {0, 1}, where the value of Vj^ [x] is determined 
by the number of particles present at x at time t. Assume that j]^ = 5o- 
Independently of each other: 

1. particles die at rate 1 without producing offspring; 

2. particles split into two at rate 1 + QjN . If split occurs at a; G Z'*, then 
one of the particles replaces the parent, while the other is sent to a 
site y chosen according to a uniform distribution on 'Qdix) (the nearest 
neighbouring sites of x\ If a newborn particle lands on an occupied 
site, its birth is suppressed; 

3. for each x,y E U^, with x — y G 9^, the values of 77^ at x and y are 
exchanged at rate (stirring). 

Just to clarify, when we say that events occur at a certain rate, we mean 
that times between events are independent exponential random variables 
with that rate. Let us also make a comment about rule 3 above. In terms of 
particles dynamics, it means that whenever exchange between sites x and y 
occurs, a particle at x (if exists) jumps to y, and at the same time a particle 
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at y (if exists) jumps to x. If one follows the motion of a typical particle, 
then, in the absence of branch events, it undergoes a symmetric random walk 
on Z'^, with jumps at rate 2dN . 

For our proofs it will be convenient to deal with the speeded-up contact 
process r}^^ : Z — t- N U {0}. This process is defined as fj^ = rj^^.t > 0. 
Clearly, r}^ = Sq. Obviously this process obeys the same rules as 77^, just all 
the events occur with rate multiplied by TV. In particular, particles die and 
split (if possible) with rates and N + 9 respectively; stirring between any 
two neighbouring sites occurs with rate A^^. 



3 Proof of Theorem 1.3 



As we have mentioned already, (11. ip and Lemma 11.41 imply that Theorem 
11.31 will follow from the lower bound for \c{N) — 1 obtained in the next 
proposition. 

Proposition 3.1 For d > 3 

N^oo ■tj/N 

In fact. Proposition 13.11 follows easily from the following crucial result. Recall 
that fj^ = Tj^iit > 0, was defined in Section [2l 

Proposition 3.2 Fix an arbitrary 9 < d. For any to > there exist Ng > 0, 
such that, for all N > Ng and t > to, 

where \ ■ \ denote the total number of particles in the process. 
Now we are ready to finish the 

Proof of Proposition [3.11 

Fix 9 < and choose Ng as in Proposition 13.21 Then, by this proposition, 
and the fact that the number of particles is an integer, we have that 
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P[\vf\=0] > 1-mf 
^ 1, 



as t ^ oo, for all > Ng. 



From this, it follows immediately that for N > Nq, {'fjf}t>o dies out in 
finite time, with probability one. The same happens with {vf}iyQ with 
probability one. 



Thus we have shown that for any 6 < -d, there exists an Nq such that for 
every N > Ng, 

P [r]f ^ 0, for all t > 0| r/^ = Sq] = 0. 

Therefore, 

inf{^ : P [7]f ^ for all t > 0| r/^ = 6o] > 0} > ^, 
and, by definition of Xc{N), the proof of Proposition 13.11 is finished. □ 



Remark 3.3 By Lemma I.4 and ( fi.ij) . we have proven Theorem \1.3\ modulo 
Proposition VJ.iA 

The rest of this section is devoted to the proof of Proposition 13. 2[ Before 
we proceed to its proof, let us first define some additional notation. 

• Denote particles by Greek letters a, /3, 7 with the convention that is 
the ancestor of a in generation 0. 

• Consider that particles change names every time an event they branch 
(die or split). 

• Let a A /3 denote the most recent common ancestor of /3 and 7. 

• In the "speeded- up" process {'t)^}t>Q, let Tq denote that time at which 
a branches (dies or slits) and let be the location of a lineage at t 
with the convention that 5" = A if the particle is not alive at t. 
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• Let (a be the indicator of the event that a spht occurred at but one 
of a's children imminently collided with another particle and died "at 
birth". 

The idea behind the proof is that in order to bound from above the total 
mass of the process we can ignore collisions between distant relatives. To this 
end, we define a sequence of times tn such that collisions between relatives 
farther related than rjv, in the process fj^ , can be ignored. 



Let 



InN , , 

rN = ^>0. (3.1) 



Let 

Z,it) = l(Ti G [t, t + r^), - e g,), 

where /3 and 7 are the children of 1 = /3o = 70 - the first particle. Zi{t) is 
the the indicator of the event that the lineage of 1 has exactly one splitting 
event in [t,t + tn), no deaths and its two offsprings (/5 and 7) are alive and 
neighbours at time t + Tpf. Set Zi = Zi{0). 

Before we proceed to the actual proof of Proposition 13.21 we will need an 
auxiliary result. Let {V/^}t>o be a continuous time, symmetric random walk 
on Z'^ jumping with rate AdN^ and starting at the origin. Let {W^}t>o be a 
continuous time Markov chain taking values in Z*^, starting from the origin, 
and evolving as follows. If Wf = x G then, with rate {Ad — 1)A^^, 
makes a jump to ?/ G Z"', whereas with probability y is chosen uniformly 
from (5d{x) \ {0}, and with probability j;^, y = —x. If W^^ = x ^ Qd, then, 
with rate AdN"^, makes a jump to y uniformly distributed in 

Note that Wf^ describes the behaviour of the difference in locations of 
two typical particles in the process fj^ in the absence of branching events. 
Such particles move around independently like symmetric random walks, 
with jumps rates 2dN'^, until they become neighbours. While they are neigh- 
bours, their behaviour is dictated by the stirring rules. 



Lemma 3.4 



E 



Jo 



= E 









W,^1{W,^ G Qd)ds 
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Proof Given V/^ G dd, define to be the time V'^ spends in 9^ before 
leaving tlie set 9^ U {0}. Tlien clearly, 

R 
i=l 

where e^s are independent random variables distributed according to 
exponential distribution with rate AdN^ and R is independent of them and 
is geometric with parameter ^f^. Clearly Ty is exponentially distributed 
with ETy = {2{2d-l)N'^)-\ 

Similarly, at every visit to S^, spends in 5^ a time = ^f^i e^, where 
R' is geometric with parameter |^5f independent exponentially 

distributed random variables with rate [Ad — 1)A^^ and independent of R'. 
Thus, is exponentially distributed with the mean 

{i^d-m^^^y = {2{2d-l)Nr'. 

Thus we may couple the processes together by setting them to be equal, 
every time they exit 9^ U {0}, and as this does not change time they spend 
in dd the result follows. □ 

Lemma 3.5 

lim NE [Zi] = d'd. 

Proof First we calculate the probability of Fi, the event that there is 
exactly one birth in I's lineage in units of time and no deaths on any of 
the branches. As the births in fj^ process occur according to a Poisson 
process with rate iV + ^, we have that 



E [z,] = p - G g,| f,] p [f,] , 
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{N + e)r^e"(^+^)^^e-^^^ • — / e-^^^'+'^'ds (3.2) 

T"7V Jo 



N + e 

2N + 9' 



-e 



-{2N+9)tn _ g-(2Ar+e)Tjv j 



where (3 and 7 are the offspring of 1 ahve at tn- 



1 fTN 

p - e F,] = — p[w + G s,] dt, (3.3) 

where W is uniform on and is the difference of positions of the two 
children of 1, right after the spht; {Wf^}t>o is a continuous time Markov 
process defined in Lemma 13.41 independent of W. 

Change the variable in the integral, use Lemma 13.41 to get that (13. 3p is 
equal to 



— / P[W + V,^ e Sd] ds (3.4) 
tn Jo 

= — y^P [W + Vne^d] -ds, 



n=0 

where is a simple symmetric random walk on Z'^ independent of W . 
Now let Ti{u) be a Poisson random variable with mean m, and define 

h{u) = P[W + e Srf] . 

Thus (13. 4p can be written as 

— / h{AdN^s)ds. (3.5) 
^iv Jo 

So, from (13. 2p . (13. 5p . and (13. ip we have that 



TJV 



EZi = E^^e-^iN+e)r^ _ ^-{2N+e)r^^ ±_ J h{AdN^s)ds (3.6) 
1 N + e r^'^^'^N 



-e 



where the last equality follows from changing the variable inside the 
integral. 
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Now let — > oo, use the Taylor expansion for the exponential and the 
monotone convergence theorem to get that 

1 

lim NEZi = — P[W + V^^s) e g,] ds. 
The times between jumps of Kr(s) are exponential with mean 1. Therefore 

^ oo 

lim NEZ,{tn) = 71 E P [^"^ ^'^l • (3.7) 

n=l 

This finishes the proof of this lemma. □ 



Proof of Proposition 13.21 

Set = E|?7j^|. From (1.5) of [6], we have 



m 



N 
Nt 



1 + 



Nt 
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N' 



ds 



2d 



Nt 



I^ds 







1 



where is twice the expected number of pairs of neighbours in rj^ at 
time s. Let A{s) be the set of particles which are alive at s in -f)^ , and 



E 



be twice the expected number of pairs of neighbours in f]^ at time s. 



Thus, with 771 



N 



771 



771 



N 
Nt 

N 



E|?7]^J, we have 



Orh^ds-—. / NV^ds. 



2d 



tN, 



Similarly to the definition of Zi, let 



where /3 and 7 are the children of a. Further, let Ca(^) be the indicator of 
the event that one of a's children created in [t, t + tn) died as a result of a 
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collision at the time of its birth Tq,. For any s > t^, let us now give a lower 
bound for 



E[Z,(s-r^)|^ 



S-TJVJ 



Note that given ^s-t^i for s > tn and a G A{s — tn), Za{s — tn) is 
stochastically less than Zi, since one of a's children created in [t,t + tn) 
may be killed as a result of a collision at the time of its birth T„. However, 
if we "return" the killed children back, then we easily get that, 
conditionally on ^s-t^, for s > tn and a G A{s — tn), 
Za{s — tn) + Ca{s — tn) is stochastically greater than Zi. Therefore we 
immediately get, that for s > tn, a G A{s — tn), 



E[Z^{s-TN)\^s-rj, 

Use this to get, for s > t^, 



^ 1 > E[Z^]-E[Cais-TN)\^ 



fN 



> E 



2E 



> 2E 



J2 l{B:-B^e<5, 



(3.9) 



E[Z^{s-tn)\^s 



J2 (E[Zi]-E[C(s-r^)|^,„.^]) 

a6(s-Tjv) 



2E [\Ais-TN)\E[Zi]]-2E 



Note now that for s < we have the same result through a similar 
calculation (with ^(0) instead of A{s — r^v)) and a slight modification of Zi. 

Using Lemma 13.51 and (13. 9p , (13. 8p gives 



(3.10) 



where Sf = NE E.^AHs-r^M U{s - tnU)- 
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On the other hand, let = s — {s — Tn)+, then 



a£A{{s-TN)+) 



(3.11) 



X— — E 

2dTt 



' N 



J2 i{b:-b:'. edd)dr 

{s-tn)+ ^^A{r) 



N 



2dT: 



■(2iV + e)r^e-(2^+^)^^E 



TV 



< CN^ I if dr. 

J {s-tn)+ 

Changing the order of integration yields 

"t rs rt /•(r+Tjv)Ai 




r:drds 



J(s-Tjv) + 




ifdsdr 



(3.12) 



Jr 



< TN I Ifdr 
'o 



where the last inequality follows from ( 13. 8p . 
From ( I33TD and (KWi we get that 

[ Efds < 2N-^tn + 2eN-^ ! mfds. 
d Jo Jo 

Thus, f l3.10p can be rewritten as 

mf < 1 + 2N-^tn + ^ (1 + 2N-^) [ mfds - U - '■ 

Jo ^ 

Clearly, from ( 13. Sp . 

Vt > 0. 



(3.13) 



rnf < 1 + ^ / m';ds, 
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Therefore, as 1 is a non-decreasing function, Gronwall's lemma gives that 

mf < e^(*-'^)mf , VrG[0,t]. 

Thus, 



By choosing A^i large enough so that for any N > Ni such that 
- e/A)e-^^^ >i^- e/2) we get. 



< 1 + 2N-Wn + {1 + 2N-^) [ mfds - -e/2) f rh^ds 

Jo Jo 

= 1 + 2N-^tn + (9 -^ + e/2 + 29N^^) [ rh^ds. 

Jo 

Again by Gronwall's lemma 

mf < (1 + 2iV-Vjv)e('-''+^/'+''^")*, Vt > 0, iV > N^. 
Since 6 < 1!} we can choose e={^-9)/2 and A^2 > A^i V to get 

mf < (1 + 2iV^ V^)et(^-'')*, VA^ > iV2, ^ > 0. 
Then choose to ^ and Ng > N2, such that 

(1 + 2A^g-V;v)e-^(''"')*° < 1, 

and hence 



□ 
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